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Abstract 

We consider a large neutral molecule with total nuclear charge Z in a model with 
self-generated classical magnetic field and where the kinetic energy of the electrons 
is treated relativistically. To ensure stability, we assume that Za < 2/tt, where a 
denotes the fine structure constant. We are interested in the ground state energy 
in the simultaneous limit Z — > oo, a — > such that k = Za is fixed. The leading 
term in the energy asymptotics is independent of k, it is given by the Thomas-Fermi 
energy of order Z 7 / 3 and it is unchanged by including the self-generated magnetic 
field. We prove the first correction term to this energy, the so-called Scott correction 
of the form S(aZ)Z 2 . The current paper extends the result of [iSSS] on the Scott 
correction for relativistic molecules to include a self-generated magnetic field. Fur- 
thermore, we show that the corresponding Scott correction function S, first identified 
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in |SSSj . is unchanged by including a magnetic field. We also prove new Lieb-Thirring 
inequalities for the relativistic kinetic energy with magnetic fields. 
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1 Introduction and results 

We consider a relativistic model of a molecule in three dimensions, where the kinetic 
energy of the electrons is modelled by the square root of the Pauli operator. The nuclei 
are fixed at positions R = (Ri, . . . , Rm) and have charges Z = (Z\, . . . , Zm), Zk > 0. Let 
Z = Ylj=k be the total nuclear charge. For simplicity we consider a neutral molecule, 
i.e. the number of electrons N is set to be equal to the total nuclear charge, N = Z. The 
particles are subject to Coulomb interaction and the electrons are dynamical. The kinetic 
energy operator of a single electron is 



T {a \A) := ^a- 2 T{A) + a~ 4 - a" 2 , (1.1) 

where a > is a parameter (fine structure constant) and T(A) is the non-relativistic kinetic 
energy operator given by 

l( — 1\ + A) (bchrodmger) . 

Here A is the magnetic vector potential generating the magnetic field B = V x A and a 
is the vector of the three Pauli matrices. Note that in the a — > limit (A) is replaced 
with ~T(A), its non-relativistic counterpart. We will treat the Pauli case (with spin-|) 
and the spinless Schrodinger case in parallel. For simplicity, we write the proofs for the 
more difficult Pauli case; the necessary modifications for the Schrodinger case are straight- 
forward and left to the reader. The operator T(A) acts on L 2 (IR 3 ,C 2 ); in the Schrodinger 
case T(A) is diagonal in the spin variables. 
The Hamiltonian of the molecule is 

H(Z, R,a,A) :=Y, {^W ~ £ T^^Tj) + £ u±~ v (1-3) 

j=l k=l ^ j<k |X? Xfc| 

where Tj a \A) acts in the Hilbert space of the j-th electron. The Hilbert space for the 
whole system is 

z 

U= /\L 2 (M 3 ,C 2 ). 
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Our units are /j 2 (me 2 ) -1 for the length, me 4 h~ 2 for the energy and mech~ l for the 
magnetic vector potential, where m is the electron mass, e is the electron charge and h~ 
is the Planck constant. In these units, the only physical parameter that appears in the 
total Hamiltonian (II. 3ft is the dimensionless fine structure constant a = e 2 (hc)~ 1 ~ It 
is known that max^ Z^a < 2/n is necessary for the stability of the system, even without 
magnetic field (A = 0). In this paper we will assume that max& Z^a < 2/n and we will 
investigate the simultaneous limit Z — > oo, a — > 0. 

For a given vector potential A, the ground state energy of the electrons is given by 

E (Z,R,a,A) := inf Spec H(Z, R, a, A). (1.4) 

The total energy with a self-generated magnetic field is obtained by adding the field energy 
(87r) _1 a~ 2 J |V x A\ 2 and minimizing over all vector potentials, 

£?o(Z,R,a) :=inf \E (Z,R,a,A) + / |V x A\ 2 }. (1.5) 
A I 87ra 2 J R3 ) 

Since the magnetic energy will always be finite, we can also assume that A G L 6 (R 3 ) (see 
Appendix of [FLL] for the existence of such a gauge), and we thus have 

V-A = } C- l ( [ A 6 ) 1/3 < [ \V®A\ 2 =[ \VxA\ 2 (1.6) 

by the Sobolev inequality, where |V <g) A\ 2 = Y^j=i \9iAj\ 2 . We will call a vector potential 
A admissible if A e L 6 (R 3 ), V<g>A G L 2 (R 3 ), and V-A = 0. Thus (O) can be reformulated 

as 

E (Z,R,a)=mf\E (Z,R,a,A) + -^ [ \V®A\ 2 }, (1.7) 
A I 8ira 2 J R3 J 

where the minimization is taken over all admissible A. 

The main question is the ground state energy in the large Z limit. The answer depends 
on whether relativistic or non-relativistic models are considered and whether magnetic 
fields are included or not. 

In the non-relativistic case without magnetic field (A = 0) the ground state energy to 
leading term is of order Z 7 ^ 3 and it is given by the Thomas- Fermi theory |LSj . The next 
order term, known as the Scott correction, is of order Z 2 and it is explicitly given by 

1 - 

*--X Z t (1-8) 

fc=i 

(the additional factor 2 is due to the spin degeneracy) and it was rigorously proved for 
atoms in [Hj [SWT] and for molecules in [IS], see also [SSJ. 

The ground state energy of the relativistic molecule without magnetic field up to sub- 
leading order (Scott correction) has been studied in [SSSj (an alternative proof for the 
special case of atoms, M = 1, was given in [FSW1] ): 
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Theorem 1.1 (Non- magnetic relativistic Scott correction |SSS] ). 
Let z = (zi,..., z M ) with z u . . . , z M > 0, J2k=i z * = and r = ( r i> • ■ • > t m) e 
with minfc^£ \ rk — re\ > r for some r Q > be given. Define Z = [Z x , . . . , Zm) = Zz and 
R = Z -1 ' 3 r. Then there exist a constant E TF (z,r) and a universal (independent of z, r 
and M) continuous, non-increasing function S2 '■ [0, 2 fir] — > R with S^O) = 1/4 swc/i i/iai 
as Z = Y^k=\ Zk ^ 00 and a — > rai/j max fc {Z fc a} < 2/n we have 

M 

E (Z, K;a,A = 0) = Z 7 / 3 E TF (z, r) + 2 ^ Z fe 2 5 2 (Z fe a) + ^(Z 2 - 1 / 30 ) . (1.9) 

fe=i 

T/ie implicit constant in the error term depends only on M and rQ. 

In the recent paper |EFS3] (see also |EFS1] and [EFS2j ) the Scott correction for a 
non-relativistic molecule in the presence of a self-generated magnetic field was proved and 
shown to be of the form S\{a 2 Z)Z 2 , i.e to depend on a through the combination a 2 Z 
(see also [Ivl llvlj IIv2j for an alternative derivation). This parameter, a 2 Z, is also the 
parameter which in non-relativistic molecules with self-generated magnetic field has to be 
small to ensure stability. Note that the physical units chosen in [ EFS 3] differ by a factor 
2 from the choice we made in this paper, in particular the non-relativistic a — > limit 
of T(A) is — in the current units, so the Thomas-Fermi energy is modified compared 
with |EFS3j . Moreover, the notation for the Scott function Si incorporates the 87r factor 
explicitly appeared in (1.8) of [EFS3J. The notations in the current paper follow the 
conventions of |SSSj . 

In the light of these previous results, it is natural to ask the following questions for 
relativistic molecules with self-generated field. 

1. (Existence of Scott term) 

Is it true that there exists a function S3 such that 

M 

E (Z, R, a) = Z 7 / 3 E TF (z, r) + 2 ]T Z 2 k S 3 (Z k a) + o(Z 2 ), (1.10) 

k=l 

in the simultaneous limit Z — > 00, a — > with k = Za fixed, for any k small? 

2. (The Scott term is non- magnetic) 

Is it true that S2 = S3, i.e. the Scott term with self-generated magnetic field is the 
same as for the non-magnetic operator (A = 0)? 

The following main theorem of this paper gives an affirmative answer to these questions: 

Theorem 1.2 (Relativistic Scott correction with self-generated field). 

Let the assumptions and notations be as in Theorem \l.l\ in particular we fix M , z and r. 

Assume furthermore that there exists k,q < 2/tc such that 

max{Zj.a} < Kq. (1-H) 
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Then the ground state energy with self-generated magnetic field is given by 

M 

E (Z, R; a) = Z 7 ^E TF (z, r) + 2 ^ Z 2 S 2 (Z k a) + o(Z 2 ) (1.12) 

k=l 

in the limit as Z — )■ oo and a — > 0. 

In contrast to the non-relativistic case |EFS3j . the Scott correction is non-magnetic in 
the relativistic case. The reason is that the prefactor (87ra 2 ) -1 in front of the magnetic 
energy is of order Z 2 in the relativistic case (since Za is bounded), i.e. it is much larger 
than in the non-relativistic case (when Za 2 was bounded). Therefore the self-generated 
magnetic field is much smaller in the relativistic case and it eventually does not influence 
the kinetic energy up to the order of the Scott term. In fact, our proof gives a somewhat 
stronger result; it proves that Theorem 11.21 also holds if the constant 8tt in (11. 7\i is replaced 
with any fixed positive finite number. 

We now comment on the new ingredients of the proof. Since the magnetic field is 
not expected to influence the final result, we can treat it perturbatively. To control this 
perturbation, our main tools are: i) a new magnetic Lieb-Thirring type inequality for the 
relativistic case; and ii) a new localization scheme for the kinetic energy operator T^ a \A). 

The magnetic Lieb-Thirring inequality for the non-relativistic Pauli operator has been 
proven in [LLS] , while the Daubechies inequality handles the relativistic case without mag- 
netic field [Dau] . Our Theorem 12.21 combines and generalizes these two classical inequalities. 
We also need a modified version of this result that allows us to include Coulomb singular- 
ities with subcritical coupling constants (Theorem 12. 3j) . We remark that magnetic fields 
have been incorporated into the Daubechies inequality even with the critical Coulomb sin- 
gularity |FLS] , but this result concerns only the Schrodinger case [FLSj where diamagnetic 
techniques are available. 

The main localization formula used in [SSSJ (Theorem 2.5) is not applicable with a 
magnetic field since it relies on the explicit formula for the relativistic heat kernel. Instead, 
we use the usual IMS formula under the square root, then apply the operator- monotonicity 
of the square root function and a useful "Pull-out" inequality (Lemma 13. ip . Constantly 
adjusting the parameter a in T^ a \A), we can show that the localization errors can be 
controlled essentially as effectively as in |SSSj despite the lack of any explicit formula. 

2 Structure of the proof 

The main steps of the proof of Theorem 11.21 follow the proof of Theorem II . 1 1 given in [SSSJ. 
To avoid unnecessary repetitions, we will sometimes explicitly refer to certain lemmas 
from [SSS] . but otherwise we keep the current paper self-contained. We will focus on the 
modifications needed in order to accomodate the self-generated magnetic field. 

We consider the number of nuclei M and the minimal distance r among the rescaled 
nuclear centers to be fixed throughout the proof and every generic constant denoted by C 
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in the sequel may depend on them. The notations 

[a] + := max{0,a} > 0, [a]_ := min{0,a} < 

stand for the positive and negative parts of a real number or a self-adjoint operator a. 
Integrals with unspecified integration domain are always considered on IR 3 . 

The upper bound in (11.121) follows from (11. 91) by choosing A = in ( 11. 51) . So we only 
need to consider the lower bound. 



2.1 Passage to the mean field Thomas- Fermi theory 

We will use the Thomas-Fermi theory for non-relativistic molecules without magnetic field 
|LS] . We will not introduce this theory here in details, we refer the reader to Section 2.7 of 
|SSS] whose notation we follow. In particular, let V£ R {x) = V TF (Z, R, x) be the Thomas- 
Fermi potential and Pz^r{ x ) — P TF (Z, R, x) the corresponding Thomas-Fermi density and 
let 

D{f,9) = \ [ I t^^dxdy, D(f):=D(fJ). 
2JJ \x-y\ 

Define the functions 

d r (x) = min {|x-r fc |}, (2.1) 

k=l,...,M 

d R (x) = min {\x- R k \} = Z- 1/3 dJZ 1/3 x). (2.2) 

fe=l,...,M 

The Thomas- Fermi potential V z T r F (x) satisfies the following bounds for all multi- indices 
nGN 3 and all x with d r (x) ^ 0: 

\%V?f(x) \ < C*mm{d T (x)- l ,d r {x)- 4 }d r {x)-W, (2.3) 

and we also have 

'<C*, for \x-r k \<r /2, k = l,2,...M, (2.4) 



K T f(x) - -^- ] 
\x-r k \ 

where the constants C* and C* depend only on r , M and m&x{Zi, Z 2 , ■ ■ ■ , Z M } (see 
Theorem 2.12 and Remark 2.14 in [SBS] ). 

We get from the correlation estimate |SSSl Theorem 2.9 (see also calculation on p. 55)] 
that if tp G H, is normalized, then 

ty, H(Z, R, a, A)ip) > Tr [T^{A) - V££(x) - CZ 3 ' 2 sg{x)\ _ 

~ D(pl%) - CsZ*l* - Cs^Z, (2.5) 

where the parameter s is chosen as 

s = Z-^ 



6 



and 

f (2s)^ 2 , d K (x) < 2s, 

g{x) = I dn(x)- 1 / 2 , 2s < d R (x) < Z^ 3 , (2.6) 

[o, Z-V3 < d n (x). 

With the above choice of s, we have sZ 8 ^ 3 = s~ 1 Z = Z 11 / 6 , so the last two terms in (12. 5p 
are negligible to the order o(Z 2 ) we are interested. 

We will estimate the error term CZ 3 / 2 sg{x) = CZ 2 ' 3 g(x) by borrowing a small 5-part 
of the kinetic energy and using the Lieb-Thirring inequality, Theorem 12.21 below (with 
h — 1 and /3 = a). With the choice of 5 = Z~ x l 2 , using a < CZ^ 1 and computing 
/ g b ' 2 < CZ- 7 ' 12 and f g 4 < CZ^ 3 , we get 



Tr 



'5T {a) (A) - CZ 3 / 2 sg(x)] > -CZ 11 / 6 - CZ 7 ' 12 ^ J |V x A\ 2 ^'\ (2.7) 



Inserting these estimates in (12.51) and using that Tr [X + Y]_ > Tr [X]_ + Tr [Y"]_ for 
self-adjoint operators X, Y, we get for any A and any normalized ip G H, 

(V, ff(Z, R, a, > Tr [(1 - zrW)T& (A) - V™] _ 

- D( P l%) - cz ll ' Q - cz 7 ' 12 ( J |v x a| 2 ) 3/4 

>Tr [{l-Z-^)T^{A)-V™]_ 

- 0(p££) " C(Z^ + Z 7 /V) - J |V x A| 2 . (2.8) 

So, using Za < C , 

W,H(Z,R,a,AW) + ^ J \VxA\ 2 >Tr[(l-Z- 1 / 2 )T^(A)-ViS_ (2.9) 

- D{pl%) - CZ** + ^—[\VxA\ 2 . 

lo7rcr J 

Clearly, the constant 16 can be replaced with any finite constant larger than 8 at the 
expense of changing C in the last line. 

2.2 Scaling 

We now introduce the usual semiclassical scaling of the Thomas- Fermi theory. The kinetic 
energy operator with the semiclassical parameter h is defined by 



[a ■ {-ihV + A)] 2 (Pauli) 
(— ihS7 + A) 2 (Schrodinger) 
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and clearly T(A) from (11.21) equals to Th=i(A). 
Define 



Za . 



k = mm 



k TTZk ' 



h = k x ' 2 Z- 1 I\ (3 = Z^an' 1 ' 2 = —h. 



(2.11) 



In particular, since Z^a < 2 /it, we have (3 < h. Note that the notation generally follows 
[SSSJ, but our definition of differs from [SSS] by a square root. 

The Thomas-Fermi potential and density satisfy the scaling relation 

Vz,n( x ) = a4v a-*z,an( ax )i Pz F r( x ) = a& P T a F ^M ax ) 
for any a > 0. In particular, 



/3 



D{ P l%) = ZV*D{p 



TF\ 
z,r I 



(2.12) 



We will perform the scaling x i— > Z 1 ^ 3 x. During the scaling we replace the vector potential 
A by 

A(x) := Z- 2 '*K l ' 2 A{Z-V*x), 
so we get for the magnetic energy in (12.91) 



\Qixoi 1 



\V®A\ 2 = Z 4 / 3 



1 



|V® Al- 



and T(A) is replaced with Th(A). Using (I2.9P and (12.121) we therefore get 

1 



(^(Z,R,«,A)^) + 



|V x A\- 



> Z^k-\1 - Z-V 2 ){Tr (^- 2 T h (A) + - /T 2 ~ l _ K z _ 1/ X r)- 

\V ®A\ 



.1/2 



+ 



16tt/3 2 /i 3 



- Z 7 l 3 D{pl F r ) - CZ U 



/6 



(2.13) 



The proof of the Scott correction is now reduced to the proof of the following semiclassical 
theorem: 

Theorem 2.1 (Scott corrected semiclassics with self-generated field). 

Suppose that A > 0. There exists a function £ : R + — > M + with £(h) — > as h — > 0, such 

that if < (3 < h, and Kmaxjzi, . . . , zm} < then 



inf {Tr y/3-*T h (A) + /3" 4 - /T 2 - W t 



TF 
z,r 



+ 



A 



/3 2 h 3 



(27ihy 



\V®A\ 

M 



-p 2 - KV^ r F (x) dxdp - 2h~ 2 Y^{zkZ) 2 S 2 {l3h- 1 



KZ k/ 



k=l 



< h~ 2 t{h). 



(2.14) 



S 



Notice that the semiclassical asymptotics up to the subleading Scott term is independent 
of the parameter A in front of the magnetic field energy. This justifies the remark that the 
specific constant 8tt in ( 11. 71) is irrelevant and can be replaced with any positive constant. 

We will first finish the proof of Theorem 11.21 using Theorem 12.11 before giving the proof 
of the semiclassical result. Using Theorem 12.11 in ( I2.13P and the choice of the parameters 
(EH}, we get 



(V,ff(Z,R, Qtj A)V> + 



57rcr 



IV x A\ 



M 



(2%h) 3 



I _ £-1/2 V *,r 



dxdp 



\ 2 S 2 [ 



(3h x ZkK 



k=l 



%)-h- 2 g{h)}-D{p TF )-CZ^ 



M 



Z 7 / 3 E TF (z, r) + 2j2 Z 2 k S 2 {Z k a) + o(Z 2 ). 



(2.15) 



k=i 



Here we used the continuity of S 2 and the facts that 

-1 . 



7^- - VTF{x) 



dxdp = C / [Vj F {x)f /2 dx < oo 



and 



z 7/3 



(27T) 



12 



p 2 -v 2 



TF I 



dxdp - Z 7 ' 3 D(p™) = Z 7 / 3 E TF (z, r) 



by standard results from Thomas-Fermi theory. This finishes the proof of Theorem 11.21 
□ 



2.3 Relativistic Lieb-Thirring inequalities with magnetic fields 

In this section we present two new Lieb-Thirring type inequalities for the relativistic kinetic 
energy with a magnetic field. The proofs are given in Section [6j 

Theorem 2.2 (Lieb-Thirring inequality for T^\A)). There exists a universal constant 
C > such that for any positive number (3 > 0, for any potential V with [V]+ G L 5//2 n 
L 4 (IR 3 ) ; and magnetic field B = V x A G L 2 (IR 3 ), we have 

Tr y$- 2 T(A) + /3" 4 - /T 2 - V{x)} _ 

> -4 J [V]f + /3 3 | [V]%+ ( | fi 2 ) 3/4 ( | [V\X) ll \ (2.16) 

Notice that Theorem 12.21 reduces to the well-known Daubechies inequality in the case 
A = [Dauj . For the Schrodinger case, the Daubechies inequality was generalized (and 
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improved to incorporate a critical Coulomb singularity) to non-zero A in |FLS] by using 
diamagnetic techniques. Theorem 12.21 is the generalization of the Daubechies inequality 
for the Pauli operator, in which case there is no diamagnetic inequality. Moreover, in 
the /3 — y limit, ( 12 . 161) converges to the magnetic Lieb-Thirring inequality for the Pauli 
operator |LLS] since 

v//3- 2 T(A) + r 4 -r 2 ^T(A), p^o. 

Theorem 12.21 does not cover the case of a Coulomb singularity. The next result shows 
that for j3 smaller than the critical value 2/n, the Coulomb singularity can be included. 
The constraint < 2/n in this theorem is the main reason why our proof of Theorem 11.21 
does not extend to the critical case, k = 2/tt. 

Theorem 2.3 (Local Lieb-Thirring inequality with a Coulomb potential). Let <p r be a real 
function satisfying supp0 r C < r}, ||</> r ||oo < 1- There exists a constant C > such 
that if (3 e (0,2 /it), then 

Tr 



(v//3- 2 T(A) + /3"4 - /T 2 - n - V)<p r 
> -C[r i/2 f I V x A\ 2 + V - 3 r 3 + tT 3/2 J [V]f + rr 3 /3 3 J [V\_ 

[v] 4 + 



+ 



IV x A 



3/4 



1/4 



(2.17) 



where rj := ^(1 - (tt/5/2) 2 ). 



For simplicity, we stated both theorems for h = 1, but T(A) = T h= i(A) can easily be 
replaced with Th(A) and the h scaling on the right hand sides can be easily computed. In 
particular, we have from ( I2.16P that 



Tr [^(3-2 Th (A)+p-*-f3- 2 -V( 



x 



> -cl hr 



[ V ]f + h -*p l[V]%+[h 



-3«3 



3/4 



[V\: 



1/4 



(2.18) 



3 Proof of Theorem 12.1 

For the upper bound in Theorem 12. II we can just take A = and apply |SSSl Theorem 1.4]. 
Notice that this will actually also provide us with a non-magnetic trial state which has the 
correct energy. 

For the lower bound we will follow the proof of the similar result in the non-magnetic 
case [SSS, Theorem 1.4] (see pages 68-75). In particular, we will use the same localizations. 
Consider a smooth partition of unity, 

e 2 _ + el = i, 
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where 0_(t) = 1 if t < 1, 0_(*) = if t > 2. Define 

$±(x) = 9±(d(x)/R), <j>±(x) = 9±(d(x)/r). (3.1) 

where := d r (x) for simplicity. At the end of the calculation we will see that the small 
parameter r can be chosen r = h 3 / 2 and the large parameter R as R = h~ l . We note that 
r was chosen differently in [SSSJ. 
For r > we define 

M 

0_(x) :=^29 rtk (x), mth9 r , k (x) := 6_(\x - r k \/r). (3.2) 

k=l 

and we note that for r sufficiently small and R sufficiently large, 9 r ^ have disjoint supports 
and these supports are in the regime where = 1. We thus have 

k 

as a partition of unity. Defining 

W V) r(x) := r _2 l{ r < d ( :r )<2 r } + R~ 2 l{R<d(x)<2R}, 

the IMS formula allows us to insert these localizations and estimate 

M 

T h (A) + f3~ 2 >J2 r,k{T h (A) - Ch 2 r~ 2 + f3- 2 )6 r , k 

k=l 

+ $_0 + (t,(1) - ch 2 w rA + /r 2 ) 

+ $ + (T,(2) - Cfr 2 ir 2 l {(i(a) < 2fl} + /T 2 ) $+. (3.3) 

Notice that all operators in the brackets in the right hand side are non-negative, since 
h 2 R~ 2 <C h 2 r~ 2 = h^ 1 <C f3~ 2 . After multiplying both sides by f3~ 2 we can take the square 
root of the inequality (13. 3p . using that the square root is operator monotone. In order to 
pull out the localization functions from under the square root, we will need the following 
general estimate: 

Lemma 3.1 (Pull-out estimate). Let I be a countable index set and let g^, % 6 I , be a 
family of non-negative smooth functions such that Yli£i9i( x ) = 1 f or ever U x £ M. 3 . Let 
A i} i e I, be a family of positive self-adjoint operators on L 2 (M 3 ,C 2 ). Then 



^giAigiy^gi^/Xigi. (3.4) 
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Proof. The proof follows from the integral representation 



y/A = (const.) J (l 



\ alt 



and from the similar "pull-up" formula for the resolvents which first appeared in [BFFGSJ, 
see also [ES.fl Proposition 6.1]. □ 

Applying the estimate (13.41) . we get from (13. 3p that 



^{3-iT h (A)+{3-*-p- 2 -KV 2 



TF 
z.r 



M I ~ 

k=l 



+ $-0+( V //3- 2 T fe (A) - Cf3-WW r , R + /3" 4 - /T 2 ~ kV^/)^- 

+ $+{ \lf3- 2 T h (A) + /3" 4 - Ch^f3- 2 R-H {d{x) < 2Ry - (3~ 2 - kV?/}*+ (3.5) 

The three terms will be considered independently in the next three subsections. The first 
one contains the contributions from the nuclei and will give the Scott term. The second 
gives the main contribution to the energy and is semiclassical. Finally the third term is a 
small error term. The main technical results, the analysis of the Scott term and the local 
semiclassical asymptotics, will be proved separately in Sections @] and [51 

3.1 The region far from the nuclei 

Let us start by considering the outer term in (13. 5p resulting from the localization procedure, 
namely 



Tr 



^ + y^T h (A) + - Ch?p-*R-n {d{x) < 2R} - /T 2 - kV™}® + ] _. (3.6) 

We introduce a dyadic partition of unity 

supp (j) , R C B(AR), supp <j> jtR C B(2 j+2 R) \ B(2 j R), 

oo 

E^ = l, \Vct> hR \<C2-=R-\ (3.7) 

j=0 



Then 



T h {A) > Y,4>3,r[Th{A) - Ch 2 2^R- 2 ]ct> hR . (3.i 
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So, 

-2rr i T\ . o-4 ^-yt2o-2 d -2 



/T^(A) + /T 4 - Ch 2 /3- 2 R- A l {d(x) < 2R} 

> E t 9 " 22 ^) - C^h^BT 3 + /3" 4 ] <j> jtR . (3.9) 



Therefore, by operator monotonicity of the square root and the pull-out estimate, Lemma I3TT1 
$ + {y / /3- 2 T,(i) + - Chif3- 2 R-n {d{x) < 2R} - f3~ 2 - 

oo , 

> E \jp- 2 T h {A) + - C/i 2 /3" 2 2-2^-2 - 0-2 - £K7}<M^- (3-10) 

3=0 

We define 

7j := 0(1 - C/i 2 /?^- 2 ^- 2 )- 1 / 4 , 

and note that < < 20 < 2/i using /3 < /i and i? > 1. We also have < /3~ 2 - 7^ 2 < 
Ch 2 2~ 2: >R~ 2 , so we can continue the estimate (using the operator monotonicity of the 
square root) as 



$ + {yJ /3-*T h (A) + (3-* - Ch 2 (3- 2 R- 2 l {d{x) < 2R} - /T 2 ~ kV£ F }$+ 

OO 

> E $ + ^«{\/77 2 ^(^) + 77 4 - 77 2 - Ch 2 2^RT 2 - kV™)®^. (3.11) 

Recall from (I2.3P that |V^ r r F (a;)| < Cd r (x)~ A and k < 2M/tt since KmaXfeZj; < 2/7r and 
maxfcZfc > 1/M. Thus \K,V^f(x)\ < C|x|~ 4 as |x| > C, where the large constant C depends 
only on M and r. Hence (V^fl < C2 _4 - 7 i?~ 4 on supp^-^. With the choice R = h^ 1 we 
can therefore absorb Vj^ in the term Ch 2 2~ 2j R~ 2 by changing the constant C. 
From the semiclassical form of the Lieb-Thirring inequality (I2.18P we get 

Tr [<!> +( j )j , R y i J 2 T h (A)+ 1 - 4 - 7 7 2 - Ch^FT 2 }^, R 

> -C{h- 3 2 3j R 3 (h 5 2~ 5j R- 5 + 1 3 h 8 2~ 8j R~ s ) + (V 2 J |V®2| 2 ) 3/4 (2- 5j ' J R- 5 /i 8 ) 1 / 4 } 

> -C{h 2 R- 2 2- 2j + j 3 h 5 2~ 5j R- 5 + 2- 2j h 5 R~ 5 + 2-^V 1 J |V <g> 2| 2 }. (3.12) 

Using the trivial estimate jj < 2h < 1 and summing up, we therefore find 

E Tr [® + ^R{\llfT h {A)+ 1 - i - 7 7 2 - Ch 2 2^R- 2 }^ R _ 
i 

> -c\h 2 R~ 2 + /i 5 iT 5 + h- 1 J |V g> 1| 2 } (3.13) 
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With the choice R = h we get 



> -C-ChT 1 / IV® A\ 2 . 



(3.14) 



So in conclusion, using the choice R— h 1 
Tr 



$ + {y / /3- 2 T,(A) + - - /T 2 - 

> -c-Chr 1 [ iv®2| 2 . 



(3.15) 



3.2 The semiclassical region 



In this section we estimate the intermediate region, the second term in (13.51) . We apply 
a multiscale analysis by localizing the intermediate regime into balls of varying radii such 
that each radius is comparable with the distance of the center of the ball to the nearest 
nucleus. We then rescale the problem in each ball to a model problem in the unit ball with 
new parameters h and /3. The model problem is analyzed in Section [5j here we state the 
scaled version of the main result: 

Theorem 3.2 (Scaled local semiclassics) . Let £, f, A > 0. Let 8 be a bounded smooth cutoff 
function supported on the ball B(£) and let V a smooth real potential on B(t). Assume 
that there is a constant C and for any multiindex n G N 3 there is a constant C n such that 

\£ H d n 6\ + \r 2 £ H d n V\ < C n , and (3f 2 £ < C'h. 



Then 



inf |Tr 



eyp^T h {A) + ^- /r 2 - v}e 
i 



+ 



Xff 



IV® A\ 



{2-nhf 



6(xY -p 2 -V(x) 



dxdp 



P 2 h 3 JB(2£) 

< CT^+i/n^-i/ii^-i/n 



(3.16) 



where C depends on X, C and on finitely many constants C n . 



Theorem 13.21 follows from the unsealed version Theorem 15.11 below with the rescaled 
variable x' = x/£ and using the parameters (3' = (3f, h! = h/(f£). □ 

The multiscale analysis requires two scaling functions, £{u) = £ u and f{u) = f u de- 
pending on u G M 3 . They express the lengthscale and the size of the potential around u, 
respectively. In our case we define 



u 



L:=^Vr 2 + d(u)\ /(«) = /„ 



mm 



{C /2 ,£- 2 }, 



14 



where we recall the definition of d(u) = d r (u) from (12. ip and that r = h 3 / 2 . The function £ u 
is essentially the distance from u to the nearest nucleus, regularized on scale r, i.e. £ u and 
d(u) are comparable if d(u) > r/3. The scaling function f 2 is the size of the Thomas-Fermi 
potential V^f near the point u. More precisely, it follows from (12.31) that 



\d n x VL F {x)\ <C n fX H , neN 3 



(3.17) 



for any x with \x — u\ < £ u and d(u) > r/3. Moreover, £(u) is a continuously differentiable 
function with ||VI||oo < 1- 

Fix a cutoff function 9 G C^°(IR 3 ), < 9 < 1, supported in the unit ball and satisfying 
/ 9 2 = 1. Define 

u {x) :=9{^) v rjM£(u) 3 / 2 , 

where J(x, u) is the Jacobian of the (invertible) map u — > (x — u)/£(u). Then Theorem 22 
from ISSI states that 



9Jx) 2 £Z 3 du = 1 



for any x G M 3 and Halloo < CJu H for any multiindex n. 

Inserting this partition of unity and reallocating the localization error, we have 



(3.18) 



f3- 2 T h (A) - CP~ 2 K 2 W r , R + /T 4 > / On P-' 2 T H (A) + /T 4 - Cf3~ 2 h 2 £~ 2 



9„, du. 



where we also used that W T , R < C£~ 2 . Since (3 < C'h and £ u > r/100 = /i 3 / 2 /100, we have 
C(3~ 2 h 2 £- 2 < Cf$- 2 hT x < /3" 4 and we can thus define 

S, := 0(1 - Ch 2 p% 2 )-^ = p[l + 0(h 2 p% 2 )} . 
Using the monotonicity of the square root and the pull-out estimate of Lemma 13.11 we get 



(3- 2 T h (A) - C(3- 2 h 2 W r , R + - /T a > / 9 U \I P~ 2 T h (A) + - ^ ~ Ch% 2 



9 U - 



(3.19) 

(Strictly speaking, the pull-out estimate was formulated for a countable partition of unity, 
but the integration over u can be approximated by a discrete sum up to arbitrary precision, 
and we neglect this technicality.) For any potential U > set 



£(A,U,9 U ) :=Tr 



$_<M„ U0- 2 T h {A) 2 -U 



ciA 



\V®A\ 



B u {2l u ) 



with a sufficiently small universal constant c%. Define the region 

Q := {u : \u\ <2R, \u - r k \ > r/3, k = 1, 2, . . . , M} 



(3.20) 
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which supports $-0+- It is easy to check that 6 I U $_0 + = for u G Q c , in particular 
£ (A, U, U ) > in this case. 

Using ( 13.191) and reallocating the field energy (this is why c\ needs to be small) we 
obtain 



Tr 



> 



where 



V+ := kVJI + Ch 2 £~ 2 . 



TF 
z.r 



A 



(3.21) 



For «e Q the localization error Ch 2 £~ 2 can be bounded by C/ 2 since /i < C£ u f u holds as 
long as c/i 2 < £ u < C/iT 1 and £ u is comparable with u. Using k < 2M/tc and (I3.17p . we 
see that 

\%V+(x)\ < C n f 2 £Z H for any x G supp(^), u G Q, n G N 3 . 
One can also easily check that 

\d n {6 u ^^ + )\<C n £-J n \ ueQ. (3.22) 

Noticing that f u £ 2 < 1 and that the condition (3 u f u £ 2 < C'h is satisfied by j3 < C'h 
(upon changing the value of C), we can now apply Theorem 13.21 to evaluate S(A, V^, 6 U ): 



(3.23) 



> 



du 



i 



[(^$_0 + )(z)] 2 [tp 2 - V +{x)\ dxdp - Ch-^fi-Wft 1 " 1 



L(2tt/i) 3 

The second term is of order /i~ 2+1 / 44 , hence negligibile, since 

/ ^l- 1/U ft 1/n < [ min{C 3 - 1/22 ,C 9+1/11 }^<C^ 1/22 <C^ 3/44 , 
is £ u Jq 

usine; that £ u ~ d(u) > r/3 = br ? ^ 2 /3. 

The double dxdp integral in the leading term is of order 

J [(^_0 + )(x)] 2 [y+(x)]f dx= J [(6 u ^<P + )(x)) 2 [kV™ + Ch% 2 ]fdx, 

which can be bounded by 

h 2 



(1 + e) / [(6 u $^ + )(x)} 2 [KV z T r F f! 2 dx + Ce 



-3/2 



\x\<CR 



5/2 



dx 



-d(x) + r. 

£ /l(^+) W ]W]>^ + C^W/< (3,4, 
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since V^f G L 5//2 (R 3 ) and i u is comparable with d(x) + r for x E supp 9 U . Optimizing for 
we see that the two error terms in (13.241) . after multiplying them with 
h~ 3 and integrating over J Q £~ 3 du, are of order h~ l ~ l l^ \ log/i|, i.e. negligible. 

In ( 13.231) we can thus replace with modulo irrelevant errors and we get 



^£(AK + A)> 



> 



du 



[($_0 + )(o;)] 



[(^$_0 + )(x)] 2 
1 



\V 2 ~ ~< r F {x) 



12 



V 



dxdp - Ch- 2+1 / u 
1 nV z T r F (x) dxdp-Ch- 2+1 '^ 



after extending the du integration to R 3 and then performing it by using (I3.18p . 
Finally, we can remove the cutoff function $_ in the last term by computing that 



(2vr/i) ; 



(1 - 



L2 



v 2 - W£ {x 



dxdp 



<Ch~ 3 j [V z T f] 5/2 l(\x\>R) 
< Ch- 3 R- 7 = Ch\ 



which is negligible. 

Thus the final result for the second term in (13. 5p is that for any A > and any admissible 
A we have 



Tr 



\^f3- 2 T h (A) - C/3- 2 h 2 W r , R + £- 4 - /T 2 - kV 2 



TF 
z.r 



+ 



A 



(3 2 h 3 



> 



(2vr/i) ; 



-p'-KV^ix) 



dxdp - Ch- 2+l / M . 



|V® A\ 2 
(3.25) 



3.3 The region near the nuclei 

Here we will consider the first term in ( 13. 5p . We consider each of the finitely many sum- 
mands individually. Without loss of generality, we may assume that = 0, so we study 



Tr 



0-{\x\/r){yl p-*T h {A) - Cf3~ 2 h\- 2 + (3~ 4 - (3~ 2 ~ ^1! )0-{\x\/r) 
A 



(3 2 h 3 



V <S> A\ 



(3.26) 



where we also borrowed a small fraction /i 1//4 of the magnetic energy. 
Define 

■= (3(1 - Ch 2 f3 2 r- 2 y 1/A > (3. 

Then we have 



< P' 2 - f3- 2 < Ch 2 r- 2 < Cti'r-'V^Xx), 



(3.27) 
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on |x| < r < ro/2. Here we used from (12.41) that for a small positive constant c (depending 
on z), we have 

rTF I 



V'*,*\x) > - C > ^ 



for Ixl < r 1. So we have 



Tr 



e-{\x\/r)(yj (3- 2 T h {A) - Cp- 2 h 2 r~ 2 + /3~ 4 - /T 2 - /^/^(MAO 



> Tr 



e.(\x\/r) V/3- 2 T,(A)+/3- 4 -/3" 



-2 ^ 



(9_(|x|/r) 



with 



// := (K + Ch'r-^izk + Cr) 



where we used the estimate 



z ' r v ' - X ~ X 



for |x| < r from ( 12.4ft . By scaling x = /i 2 /i 1 y, this becomes 
Tr 



0_( M/r) ( V //3- 2 ^(^) + /3" 4 - r 2 - )e_(H/r) 



1 



with 



Notice that 



therefore 



0-(\y\/K) ( x /5- 2 ^=i(A) + 5-4 - 5T 2 - — )0-(\y\/H) 



H / |V® A{y)\ 2 dy = / | V <g> A{x)\ 2 dx, 



(3 2 h 3 



\V ® A(x)\ 2 dx = ^{k j \V®A{y)\ 2 dy]. j 



with 



A = ^ >>L 



So with the choice r = h 3 ^ 2 we have (using (3 < h and KZk < C), 

H 5 /A < Ch 1/4 , 



(3.28) 



(3.29) 



(3.30) 



(3.31) 
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where the constant depends on A. With this choice we also have 

a = /3/rV = (5h-\K + Ch 2 r' l )(z k + CV)(1 - Ch 2 p 2 r~ 2 )- 1/4 = f3h- 1 nz k + 0{h 1/2 ). 

So using Lemma H~T1 below and the continuity of S 2 , we get, after rescaling to the original 
coordinates that 



lim inf h 2 { Tr 



6-(\x\/r) yP~ 2 T h (A) - C(3- 2 h 2 r-z + /3" 4 - (3~ 2 - MM AO 



A 



IV® A\ 



/W J 1 1 (2vr/i) 3 

2(^ fc ) 2 S 2 (/3/i- 1 ^ fc ) I > 0. 



9 2 _{\x\/r) 



12 



2 

v - 



s J - 



dxdp 



(3.32) 



We can replace /VM with kVJ t (x) in the semiclassical formula, the error is of order 





r n ] 


5/2 


J Bl{\x\/r){ 








-\x\. 





X 



< + ^3/2 < ^5/4 



which is negligible, where we used (j2.4p and that /x = KZk + 0{h 1 / 2 ). After this replacement, 
we can sum up (I3.32p for each k to obtain the final result of this section: 



M 



lim inf h 2 > < Tr 



T fe 



'r,k 



k=l 



A 



P 2 h 3 J 1 1 (2tt/i) 3 
- 2fe) 2 S , 2 (/3/r 1 £z fe ) j> > 0. 



z,r 



dxdp 



(3.33) 



Combining the estimates A3. 15j) . (I3.25j) and (I3.33P on the three terms in (13. 5 j) and 
recalling (13. 2p . we immediately obtain (I2.14p . This completes the proof of Theorem 12.11 
□ 



4 The Scott region 



In this section we fix a non-negative cutoff function <fi : M 3 — > [0, 1] with support on the 
unit ball -B(l) and such that (f> = 1 on 5(1/2), the ball of radius 1/2. Set 4> r {x) : = (fi(x/r) 
for any r > 0. Define, for R, A > and a G (0, 2/ir), 



£ w (A)=Tr [<f> R (Va- 2 T h =i(A) + or 4 - a~ 2 - -^)cj> R \ _ + A / |V<M 



(27T) 5 



1 

I .r I 



dxdp 



(4.1) 
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and 



E(R,a,A) = inf£ RiaA {A). 



(4.2) 



Clearly, E(R,a,A) < £ r ^k{A = 0), and we know from |SSS|. Lemma 4.3] that 
£r,o,a(A = 0) tends to the non-magnetic, relativistic Scott term 2S 2 (a) (the factor 2 
being due to the spin degrees of freedom). 

Lemma 4.1. Fix £ (0,2/tv). We take the limits R, A — > 00 in such a way that R 5 / A — > 
0. Then we have 



lim E(R,a,A) =2S 2 (a), 



(4.3) 



where the limit is uniform in a < ojq. 



Proof. As mentioned above, the upper bound follows by taking A = and using [SSSl 
Lemma 4.3]. We proceed to give the lower bound. 

Step 1: A priori bound on the field energy. Theorem 12.31 with V = yields 



Tr 



cf )R ya^T h=1 {A) + 



a-*-a- 2 -l/\x\)<p R 



> -C 



IV ® A\ 2 + R' 



(4.4) 



with a constant C that only depends on - — a > 0, i.e. on the distance of a from its 
critical value -. Notice that the Weyl term also satisfies a similar bound, 



<t> R {x) 



1 2 K 

2 P ~ R 



dxdp 



C^' 2 R l '\ 



(4.5) 



for some constant only depending on (p. 

Inserting these bounds in £ Ria ,A we get for any A with £ Rt a,A(A) < £ Rt a^{A = 0) that 



(A - C) J |V ® A| 2 - C(i? 3 + i? 1 / 2 ) < ^, Q , A (A = 0). 



(4.6) 



We know from [SSSJ that £R >a) A(A = 0) tends to 252(a) and 5*2 is bounded by 1/4. In 
particular, the right hand side of (14.61) is bounded by some constant C for large values of 
R. So we get for all A, R sufficiently large that 



|V<8> A\ 2 < CR 3 /A. 



(4.7) 



Step 2: Localization of the vector potential. We start by localizing the vector 
potential A. Suppose that A satisfies that £R )a ,\(A) < £ Rt a,A(A = 0). We may add a 
constant to A without changing £ R . a ,h{A). So we will assume that 



Adx = 0. 



B(2R) 
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Let X11X2 £ C°°(M ) be a partition of unity satisfying 

X?+X2 = l, Xi = lonB(l), suppxiC J B(2). (4.8) 

Define Xj.-RtX) = Xj{ x /R)i 3 = 1)2. Let also xi £ Cq°(B(2)) with %i = 1 on suppxi- We 
define 

M x ) = Xi( x / R )A( X )- 
With this notation we get from the IMS-formula (and since X1X1 = Xi) 



a 



~ 2 T h=1 (A) + a" 4 > X i, R [cx- 2 T h=l {A) - Ca~ 2 R~ 2 + a^]xi,R 

+ X2,n[a- 2 T h=1 (A) - Ca 2 R~ 2 + a" 4 ]^- (4.9) 

Using the operator monotonicity of the square root, the pull-out estimate of Lemma 13.11 
and that §bX\ r = &R, 4>rX2 r = 0, we therefore have 



a- 2 T h=1 (A) + a-^ R > <j> R d a^T h=l {A) - Car*R-* + a~\ 



'R 



> <p R ^-*T h=l {A) + 7-V fl , (4.10) 

where 

7 = a(l - Ca 2 R~ 2 )~ 1/A > a. (4.11) 
Step 3: Removing the magnetic field. To continue the lower bound, we estimate 

T h= i(A) > -(1 - 2e)A + e( - A - e- 2 i 2 ) (4.12) 

with some e G (0, 1/2) to be determined later. By the CLR-inequality 

-A - C 2 A 2 > 

if 

Cclr J {e~ 2 A 2 f 2 < 1, (4.13) 

where Cclr is an explicit constant in the CLR inequality. By the Cauchy-Schwarz and 
Sobolev inequalities, and using that A is supported on B(2R), we obtain 

J(e- 2 A 2 f 2 = e- 3 ( J l) V2 ( J A 6 ) 17 ' <Ce- 3 R 3 / 2 ( J |V®i| 2 ) 3/2 . (4.14) 
We can continue the estimates using the Poincare inequality (since J b ^ 2 r) Adx = 0). 
/ |V ® A\ 2 < [ |V <g> A\ 2 + \Vxi{-/R)\ 2 A 2 < [ |V ® A\ 2 + CR~ 2 [ A 2 

J J J Jb(2R) 



<C / IV® A\ 
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So we may replace A by A in (I4.14p . 

Upon inserting this estimate in (I4.14p and using (14. 7p . we see that the condition ( 14. 13[) 
is satisfied if we take 



(4.15) 



with a sufficiently large constant fi. Clearly e G (0,1/2) can be achieved in the limit 
considered in (14. 3p . 

With the choice of e from (14.151) and using (14.121) and the operator monotonicity of the 
square root, we have 



£ R ,a,A(A) > Tr [<^,(vV 2 (i-2e)(-A)4 . 
2 



-4 



a 



l/\x\ 





"1 2 


1 " 




2 P ~ 


\x\ 



dxdp. 



(4.16) 



We perform the scaling y — (1 — 2e) 1 l 2 x in order to absorb the factor (1 — 2e). With the 
new parameter 

R = (1 - 2er 1 / 2 i?, 



we get 



Tr 



fi ( v / 7- 2 (l-2e)(-A) + 7 - 



4 - a" 2 



l/|x|)0 



Tr 



g f > /r- 2 (-A) + T 4 -7- a -( 



2 -2> 

a -7 ) 



Vl^2^ 



(4.17) 



We use the Lieb-Thirring inequality Theorem 12.21 (in this case Theorem 12.21 is the usual 
Daubechies inequality) to control the difference (a~ 2 — 7~ 2 ). For this we will use a small 
5-part of the kinetic energy (in the end we will make the optimal choice 5 = R^ 1 ). Since 

< a' 2 - 7 ~ 2 < C7T 2 , 



we get 



Tr 



<^WV7- 2 (-A) + 7 - 4 - 7~ 2 ) - («~ 2 - r 2 ))^ 
> -c 



r 3/2 (a -2_ 7 -2 )5 /2 +7 3 5 -3 (a _ _ fy 



-2\4 



a -2 _ -2^5/2 , „:-!.s-o, 

'{[*[<*} 

> -C(<T 3/2 ir 2 + <T 3 iT 5 ) 

With the choice S = R^ 1 this term is estimated as CR^ 1 ^ 2 . 

For the main term, containing (1 — <5)-part of the kinetic energy and the Coulomb 



(4.18) 



potential, by scaling x 



1-6 



y, we have 



1 - 6)Tr k(v / 7- 2 (-A)+ 7 - 4 - 7~ 2 ~ ^ ~^~T~T~~~ )$r 



\x\ 



1 - 2e. 



-Tr 



^(V7- 2 (-A) + 7- 4 -7" 2 



1 

I x I 



(4.19) 
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with R = R^Er and 7 = 7^^ 



<5 • 



Notice that the classical (Weyl) terms satisfy ( 14. 5 ft and therefore 

1 - 2e 







"1 2 


1 " 




//< 




f 


\x\ 






R 1 ' 2 - 


1 - 


2e R 1 ' 2 


1 - 


-5 





dxdp — 



4>\{x) 



1-5 

<CR 1/2 [e + 5} = o(l) 



1 2 1 
2 P ~- 



x 



dxdp 



(4.20) 



using the choice 5 = R 1 and by the choice of e and since i? 1//2 e — > under the limit taken 
in g3D. 

So we can insert the above estimates into (I4.16P to get 



1 -2e 



E R , a AA) > ——£ R ~ A (A = 0) - CR l ' 2 [e + 5]- C{5^' 2 R~ 2 + 8^RT b 



1-5 



(4.21) 



Since we know from |SSS] that £rx^{A = 0) — » 262(0;) this finishes the proof using the 
previously mentioned choices 5 = i? -1 and e from (I4.15p . □ 



5 Local semiclassics 

Theorem 5.1. Let 6 be a bounded cutoff function supported on the unit ball B(l) and V 
a smooth potential on B(l). Let A > be fixed. Assume that there is a constant C and 
for any n6N 3 there is a constant C n such that 

\d n V\ < C n , and (3 < C'h. 

Then 

inf {Tr [O{y/0-*T h (A) + p-*- /T 2 - V}6] _ + J | V ® A\ 2 ) 

< Chr 2+1 ' n , (5.1) 
where C depends on A and on finitely many constants C n . 

Proof. The upper bound follows by choosing A = and applying [SSSj, Theorem 4.1]. 
Also, using Theorem 12.21 it suffices to prove the estimate when h is sufficiently small, say 
h < 10" 11 . 

Let I be a parameter satisfying h < I < 1/10. At the end of the proof we will choose 
i = h 1 ' 11 . Let {^} i6 %i be a collection of smooth functions satisfying 

4>le = 1, supp <f> jJL c B je {£), I V ^'| 2 ^ Ci ~ 2 > 

j 3 



(2nhy 



9{x) 



-p 2 - V(x) 



dxdp 
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where B x (r) is the ball of radius r centered at x. Then, by the IMS-formula 

T h (A) > J2he(Th(A) - Ch 2 r 2 )<P 3/ . 
j 

So by the pull-out estimate of Lemma [3.11 and operator monotonicity of the square root, 



Tr \eyp^T h (A) + {3-±-{3-*-V}6 
> ^Tr [0<P hi yi3^T h {A)-C(3 



(5.2) 



- 2 h 2 £~ 2 + /3" 4 - (3~ 2 -vu ie e 



Also, with some universal constant cq 

j:B je (£)nB(l)^d) ^ B M 2£ ) 



V / |V® A\ 2 < c / |V® A\ 2 , 

JB{2) 



SO 



Tr \ey[3- 2 T h (A)+(3- 4 -p- 2 -V}6 



+ 



A 



IV® A\ 



2 h 3 j B{2} 

> { Tr [0<l>jA VP- 2 T h (A) - Cf3- 2 h 2 i- 2 + (3-* - fl~ 2 - V}4> 3/ 9 
j 

A 



+ 



c (3 2 h 3 



(5.3) 



B ]t (2t) 



We may consider each summand independently and we can focus only on those that give 
negative contribution. For simplicity of notation, we take j = 0. Choose a new, smooth 
partition of unity x\ + x\ = 1> with 

suppxi C B(2£), xi = 1 on supp0 O /, |Vxi| 2 + |Vx 2 | 2 < CT 2 . 

By a constant shift in A and gauge invariance we may assume that § B , 2 (\ ^4 = 0. 
Choose A = xA, with x satisfying the same conditions as Xi an d XXi — Xi- Define 
^ M = / i 2 (|V X i| 2 + |V X2 | 2 ). 

Then, by IMS and the pull-out estimate again 

4>^f3- 2 T h {A) - Cf3- 2 h 2 £- 2 + p-% >e 



M f3- 2 (XiT h (A) X i + X2T h (A) X 2 - W h , e ) - C/3~ 2 h 2 e~ 2 + 



> <j> jtt y/ /3- 2 T h (A) - C'l3- 2 h 2 £~ 2 + 
with a different constant C, where we used that Xi^j/ 



(5.4) 
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Using the Poincare inequality, we also have 

I |v®l| 2 = / \V®A\ 2 < [ x 2 |V® A\ 2 + 2\Vx\ 2 A 2 

JR 3 JB(2l) Jb(2£) 

12 



< d / |V® A\ 2 (5.5) 

for some universal constant C\. 

So for each j, it suffices to consider a semiclassical lower bound to 



inf < Tr 

A 



i jtt {y/l3- 2 T h (A) - cp-we- 2 + p-* - /r 2 - v}^ 

+ \V®A\ 2 }. (5.G) 



c /3 2 /i 3 



B j£ (2£) 



where cq is a given fixed constant, and the infimum is taken over all vector fields A with 
support contained in Bji(2£). 

First we get a crude upper bound on the field energy. Clearly 



y/p- 2 T h (A) - C(3- 2 h 2 £- 2 + /3" 4 - f3~ 2 > ^p- 2 T h {A) + /3~ 4 - ft' 2 - Ch 2 £~ 2 
with (3 defined by 

p-* = /T 4 - C(3- 2 h 2 £~ 2 . 

We can apply the Lieb-Thirring inequality Theorem 12.21 for the first line of (I5.6P with 
the bounded potential V(x) + Ch 2 £~ 2 l(x e supp <pjj) to obtain a lower bound of order 
—Ch~ 3 £ 3 . This implies that 

B 2 := [ \V®A\ 2 < C[3 2 £ 3 , (5.7) 

JB(2£) 

with a large constant C whenever the vector field A gives a non-positive energy. 
We now estimate, for any e 6 (0, 1), 

T h (A) > -(1 - 2e)/i 2 A + e(-h 2 A - e~ 2 A 2 ) 

>-(l-2e)h 2 A-Ch- 3 e-H 1/2 B 5 . (5.8) 

Here the last inequality follows from the Lieb-Thirring, Holder and Sobolev inequalities 
recalling that A is supported in B^{21): 

-h 2 A - e- 2 A 2 > -Ch- 3 e- 5 [ A 5 > -CTrV 5 (I 1 V'Y / A^ '" " 



B(2£) 



> -Ch~ 3 e- 5 £ 1/2 ( J \V®A 



5/2 
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Define 

7 - 4 : = /T 4 - C(5~ 2 hH- 2 - C(3- 2 h-h-H l ' 2 B\ h = hVT^Ye. (5.9) 
We will in the end make the (optimal) choice 

e = h- 3 ' 6 £ 1 ' w B. (5.10) 
Using (15.71) . this choice will ensure that 

e < l, h-h-H^B 5 = h-Vn i ' 10 B < 1, 
so 7 is well defined. Using operator monotonicity of the square root we get from (I5.8P that 
^f3-*T h (A) + (3-±-(3- 2 -V{x) (5.11) 



> V7" 2 (-^ 2 A) + 7" 4 - P ~ V(x) 



> \jl- 2 {-h 2 &) + 7" 4 - 7~ 2 - {V(x) + Ch- 3 e~ 4 £ 1/2 B 5 + Ch 2 £~ 2 ). 
Using [SSSl Theorem 4.1] we therefore have 

A 



Tr 



> 



2 



+ 



c (3 2 h 3 



V <g) A\ 



B(2£) 



, ^x) 2 2 ,(x) vW+F-r 2 -^) ^-cw^r 9/5 

(ZTTAij' 3 7 J ' L J - 

C(/i/£)- 3 (e + h- z e-H l ' 2 B h ) + /T 2 /^ 2 . 



(5.12) 

The leading semiclassical term is of order (h/£) 3 . The e term in the second line comes from 
adjusting h to h in the main term and we have absorbed the term C(h/£)~ 3 h 2 £~ 2 (from the 
last line of (15.1 ip ) in (h/£)~ 9 ^ 5 . In the leading term we can replace v^V+F _ p-2 
with ^p 2 at the expense of a (3 2 {h/£)~ 3 error. 
By the choice of e above, the last line is 

-Ch- 3 £ 3l ' w h- 3 ">B + f3- 2 h- 3 B 2 . 

Clearly, 

Ch- 3 £ 31 / 10 h- 3 / 5 B < (3- 2 h- 3 B 2 + C 2 h- 3 £ 31 '^ 2 h-^ < (3~ 2 h- 3 B 2 + C 2 h- 3 {£ 3l ' 5 h^), 
using the bound (3 < C'h. So we get 

A 



Tr 



*o4V 'P- 2 Th{A) + /3" 4 - /T 2 - U(x)}0 o ,^ 



+ 



> 



^(a;) 2 ^ 2 ^) y-p 2 - V(z)J dxdp - C(/i/£r 9/5 



(2tt/i) 3 

-C(V^)~ 3 (/3 2 + ^ 16/5 /i 4/5 )- 



(5.13) 
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With the choice £ = h 1 ' 11 , we have 

(h/£)- g/5 + (V^r 3 (^ 16/5 fr 4/5 ) = 2(h/£)- 3 h 12 / u 

and the error term from (3 2 is negligible since (3 < C'h. 

Similar bound holds for any j. We proceed to insert (15.13!) for each j in (15. 3ft and get 



Tr 



6{y/^T h (A)+^- /T 2 - V{x)}( 



A 



c fi 2 h 3 J B{2) 



IV® A\ 



> 



9{x) 



v ~V 2 -V{x)\ dxdp-Ch- 3+12/n . 



(5.14) 



Here we used that the summation in (15.31) can be restricted to those j, where the ball 
Bj£(£) has non-empty intersection with B(l). By a volume argument there are of order of 
magnitude £~ 3 such balls. 

□ 



6 Proof of the relativistic Lieb-Thirring inequalities 

Proof of Theorem \2. 6 A By scaling it suffices to prove that there exists a constant C > 
such that for all m > 0, 



Tr 



\/T(A) + m 2 -m- V(x) 



> -C{m^ J [V]T + J [V]% + (/ |V x A\ 2 f\ j 



[V\: 



(6.1) 



The basic idea is to consider the spectral subspaces on which T = T(A) < 10m 2 and its 
complement, T > 10m 2 , separately. For any T > and m > we have the following 
simple arithmetic inequalities: 



VtT 



m > 



c T 



ni 



if T < 10m 2 



VT + m 2 - m > -Vf, ifT>10m 2 
3 



(6.2) 



where c is a universal constant (actually c = ^(vTT — 1) will do). In the first regime 
we can use the non-relativistic magnetic Lieb-Thirring inequality. In the second regime we 
will use the BKS inequality |BKS] stating that 



Tr (P - Q)_ > -Tr |-(P 2 - Q 2 )_ 



1/2 



(6.3) 



for any positive operators P, Q. In this way we can turn the problem in the second regime 
into a Lieb-Thirring type estimate on the half moments of the negative eigenvalues of the 
Pauli operator. 
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For the detailed proof, we can clearly assume that V > 0. We split the potential as 

V = Vl + V 2 , Vy{x) := V(x) ■ l{V(x)<m/2}, V 2 (x) := V(x) ■ l{V(x)>m/2}- 

Since 



Tr 



\JT + m 2 — m — V\{x) — V 2 (x) 



> -Tr 
~ 2 



\/T + m? — m — 2Vi(x) 



+ ^Tr 
2 



VT + m 2 — m — 2V2(x) 



for the proof of (16.11) it suffices to consider separately the cases 

• V(x) < m for all x; 

• V(x) > m whenever V(x) ^ 0. 

The first case will be applied to V being 2Vi, while the second to V being 2V2. 
In the first case, V < m, consider the projections 



:= 1 



< : ~ J -{T<10m 2 }, P> '■— l{T>10m 2 }) 



and estimate 



V = (P< + P>)V(P < + P>) < 2P < VP < + 2PyVP>. 
Since V < m, we have 



VT + m 2 -m-V > P< (Vf + m 2 - m - 2V\P< + P> (y/T + m 2 - m - 2V^jP> 

> P < (m- 1 c T -2V^)P < , (6.4) 



where we used the spectral theorem and the elementary inequality (16. 21) . 
Therefore, 



Tr 



VT + m 2 — m — V(x) > m x Tr [coT — 2mV]_ 

>-C{m^j[V]r+(jBf\J[V] 



(6.5) 



1/4 



where the last inequality follows by the Lieb-Thirring inequality for the Pauli operator. 

We now consider the case where V(x) > m whenever V(x) 7^ 0. Notice, that in the 
special case m = this condition is automatically satisfied. Here we first use the BKS 
inequality (16. 3 j) to get 



Tr 



VT + m 2 — m — V(o 



> -Tr - 



> -Tr 



T + m 2 - (m + V) 2 ) 
T-3V 2 



1/2 



1/2 



(6.6) 
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Here we estimated mV < V 2 by the assumption on V to get the last inequality. 

We now use the "running energy scale" method from |LLSj . For a self-adjoint operator 
H let N{H) denote the dimension of the negative spectral subspace. Let A G [0, 1] be a 
real parameter chosen at the end. Using that T > and T > (p + A) 2 — | V x A\, we obtain 



Tr 



T-3V 2 



1/2 



M{T-3V 2 + e)^ 
o V e 



/>oo 

< / A/"(AT - 3V 2 + e) 

/>oo 

< / M{\{p + Af -\\V x A| -3V 2 + e) 



(6.7) 



We now apply the CLR-estimate to get 



Tr 



T-3V Z 



1/2 



< c 



de 



[|V x A| + 3A~V 2 - A^e] 3 / 2 ^^ 



<C'{V\ J |V x A\ 2 + \- 2 1 V 4 ]. 



(6.8) 



Setting S = |V x A| for simplicity, if J B 2 < jV A , we choose A = 1 and get a total 
estimate of size J V A . If J B 2 < J V A we choose A = (J U 4 / J B 2 ) 1 / 2 and get an estimate 

/ \ 3/4 / \ 1/4 

of size f / B 2 ) (J V A \ . Thus the choice A = min {l, (/ V A / J B 2 ) 1 ' 2 } will do the job 
in both cases. This finishes the proof of (16.11) and therefore of Theorem 12.21 □ 



Proof of Theorem \2.3\ . A potential with Coulomb singularity is not allowed in Theorem [221 
We will need to use Kato's inequality to control the Coulomb singularity directly; the 
remaining part of the potential will be treated as in the proof of Theorem 12.21 This will be 
done in Lemma 16.11 However, this estimate does not have the expected behavior for small 
values of (3, where one should be close to the non-relativistic situation. So Lemma [6. II will 
be used only if (3 is separated away from zero, say (3 > 1/20. When (3 G (0, 1/20), we can 
estimate ((3~ 2 T + /3" 4 ) 1 / 2 - f3~ 2 by (const. )T effectively, and we use a result from |ES2] to 
"pull the Coulomb tooth" . This is the content of Lemma 16.21 below. 

Lemma 6.1 (Stability up to the critical coupling). Let r > r for some given r > 0. Let 
4> r be a real function satisfying supp0 r C {\x\ < r}, \\4> r \\oo < 1- There exists a constant 
C > 0, depending only on r , such that if (3 G (0, 2/7r), then 



Tr 



<p r ^Jf3- 2 T(A) + (3-±-f3- 2 - — -V 
V \x\ 



v x a\ 2 + /rV 3 ^ 3 + rr 3/2 



W)¥+ti-*0* im 



3/P3 



+ 



|VxA| 2 ) 



3/4 



/3 -l/2 r 3/4 + 



[VI 



1/4- 



(6.9) 



where r\ 



4(1 



(7T/3/2) 2 ). 
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Proof of Lemma \6. 1\ Without loss of generality we can assume that V > 0. We estimate 
<p r yp-*T(A) + f3~± - /T 2 - jl - \/]0 r (6.10) 

> 4> r [(1 - rj) VP~ 2 T(A) ~ T^j] 0r + r h V^TpyT^ 1 - /T 2 l {N <r} - 0r- 

In the first term we use the Kato inequality (2/tt)/\x\ < | (— iV+A) | and the BKS inequality 
U53D to get 



TV 



1/2 



1/2 



> Tr [0 r ((l - v ) V(3- 2 T(A) - (7r/2)|(-iV + A)|)0 r 

> Tr [(1 - riWp-*T{A)- (7r/2)|(-iV + A)| 

> -Tr (- [(l- r? ) 2 /3- 2 T(A)-(7r/2) 2 |(-zV + A)| 2 ] J 

> -/3^Tr ( - [((1 - r?) 2 - (/3tt/2) 2 ) (-iV + A) 2 - (1 - V ) 2 \ V X A|] 

> -/r lr Tr ( - [8?7(-zV + A) 2 - |V X A|]_J 

> -Cf3~ l il~* 12 j |V x A| 2 , (6.11) 

where we also used the trivial lower bound for the Pauli operator, T(A) > (— zV + A) 2 — 
|V x A|, in the fourth line and the special choice of rj in the fifth line. The last inequality 
in (16. lip is the non-relativistic Lieb-Thirring inequality for half moments of the negative 
eigenvalues. Note that the bound on this term is consistent with (16.91) . 

We proceed to estimate the second term in (I6.10p using the Lieb-Thirring inequality of 
Theorem 12.21 



Tr 



4>r(vV (3~ 2 T(A) + (3-* - /T 2 l {N < r} - V)ct> r 



> _cJ / (/r 2 i {N <n + vf 2 + /3V 4 / (/r 2 i {N < r} + v) 



IV x A 



2 ) 3/ V 1 (/(r 2 i { N< r >+^) 4 ) 



1/4 



Elementary calculations show that this is also consistent with (16. 9j) . 



(6.12) 



□ 



Finally, the expected behaviour for small values of f3 can be obtained by the following 
modified version of Lemma 16.11 
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Lemma 6.2. Suppose j3 G (0, 1/20), andr\ > 0. Then there exists a constant C depending 
only on r\ such that for any r < r% we have 



Tr 



^{3-*T(A)+(3-* -(3- 2 - — - l {N < r} - V 



> -C{l + J [V]f + J [V] 4 + + J |V x A] 2 }. (6.13) 



Proof of Lemma 1 6. 2[ Assuming again V > 0, we first split the energy as follows 

1 



Tr 



s/P~ 2 T{A) + f3~* - 13- 



\x 



M\x\<r} - V 



> -Tr 
~ 2 

1. 



+ -Tr 
2 



V/3" 2 T(A) + /3" 4 - /T 2 - 2V 
y/p- 2 T(A) + 0=* - /T 2 ~ ^iM\ x \<r} 



(6.14) 



The desired estimate for the term with V follows from Theorem l2.2[ so it suffices to consider 
the second term with the Coulomb potential. Following the proof of Theorem 12.21 we split 
this as 



Tr 



y/P~ 2 T(A) + /3" 4 - /T 2 - — 1 



x 



{\x\<r} 



> -Tr 

~ 2 

+ ^Tr 
2 



y/p- 2 T(A) + W A ~ r 2 - T^-{*P<\*\<r} 



V/3- 2 T(A)+/5" 4 - (3~ 2 - T-rl^m 



(6.15) 



The first term on the right, where the potential is bounded by (3 2 , we can estimate 
similarly to ( I6.4p as 



Tr 



^f3~ 2 T{A) + /3" 4 - (3~ 2 - — 1 



kt: 



{4/3 2 <M<r} 



> Tr 



8 

CqT(A) — — l{|x|<r} 



(6.16) 



This is a Pauli operator with a Coulomb singularity and is known to be bounded from 
below by a constant depending only on the upper bound on r, see the (proof of) Lemma 
2.1 in |ES2] with the choice of Z = Sc^ 1 (see also |EFS3t Equation (4.8)] with a special 
choice of the constants). 
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For the second term in (16.151) we use the BKS inequality ( 16. 3p and estimate 



Tr 



(6.17) 



> -Tr 
= -Tr 

> -Tr 

> -Tr 



'p~ 2 T(A) + /r 4 - (/r 2 + Ai {N < 4/32} ) 2 



1/2 



^T{A)-^-l { ^ } -- 



16 



1/2 



/t 2 t(a) - 



48 



-{|x|<4/32} 



1/2 



400{(-iV + A) 2 - |V x A\) 



48 

|x| 2 J - 



1/2 



where we used /3 < 1/20, T(A) > and that T(A) > (-z'V + A) 2 - |V x A|. We now use 
the Hardy inequality, (4|x|)~ 2 < (— zV + A) 2 to continue this estimate as 



> -Tr - 



208(-zV + A) 2 - 400|V x A\ 



1/2 



> -C / |V x A| : 



(6.18) 



by the non-relativistic Lieb-Thirring inequality for the half moments. This finishes the 
proof of Lemma 16.21 □ 



As explained previously, the results of Lemma 16.21 and Lemma 16.11 combine to imply 
Theorem 12.31 Therefore the proof of Theorem 12.31 is finished. □ 



References 

[BFFGS] L. Bugliaro, C. Fefferman, J. Frdhlich, G. M. Graf and J. Stubbe: A Lieb- 
Thirring bound for a magnetic Pauli Hamiltonian, Commun. Math. Phys. 187, 
567-582 (1997) 

[BKS] M.S. Birman, L.S. Koplienko and M.Z. Solomyak: Estimates for the spectrum 
of the difference between fractional powers of two self adjoint operators. Soviet 
Mathematics, 19(3) 1-6, (1975) 

[Dau] I. Daubechies: An uncertainty principle for fermions with generalized kinetic 
energy. Commun. Math. Phys 90 (1983) no. 4, 511-520. 

[EFS1] L. Erdos, S. Fournais, J. P. Solovej: Stability and semiclassics in self- generated 
fields. Preprint: arxiv.org/1105.0506 

[EFS2] L. Erdos, S. Fournais, J. P. Solovej: Second order semiclassics with self- 
generated magnetic fields. Preprint: arxiv.org/1105.0512. To appear in Ann. 
Henri Poincare. 



32 



[EFS3] L. Erdos, S. Fournais, J. P. Solovej: Scott correction for large atoms and molecules 
in a self- generated magnetic field. Preprint: arxiv.org/1105.0521. To appear in 
Commun. Math. Phys. 

[ESI] L. Erdos, J. P. Solovej: Uniform Lieb- Thirring inequality for the three dimen- 
sional Pauli operator with a strong non-homogeneous magnetic field. Ann. Henri 
Poincare 5 (2004), no. 4, 671-741. 

[ES2] L. Erdos, J. P. Solovej, Ground state energy of large atoms in a self- generated 
magnetic field. Commun. Math. Phys. 294, No. 1, 229-249 (2009) 

[FLL] J. Frohlich, E. H. Lieb, and M. Loss: Stability of Coulomb systems with magnetic 
fields. I. The one-electron atom. Commun. Math. Phys. 104 251-270 (1986) 

[FLS] R. L. Frank, E. H. Lieb, R. Seiringer: Hardy-Lieb-Thirring inequalities for frac- 
tional Schrodinger operators. Journal AMS 21, No. 4, 925-950 (2008). 

[FSW1] R. L. Frank, H. Siedentop, S. Warzel: The ground state energy of heavy atoms: 
relativistic lowering of the leading energy correction. Commun. Math. Phys. 278 
no. 2, 549-566 (2008) 

[H] W. Hughes: An atomic energy bound that gives Scott's correction, Adv. Math. 

79, 213-270 (1990). 

[Iv] V.I. Ivrii: Local trace asymptotics in the self-generated magnetic field, Preprint. 

larXiv:1108.4T08l 

[Ivl] V.I. Ivrii: Global trace asymptotics in the self-generated magnetic field in the 
case of Coulomb-like singularities Preprint. larXiv:1112.2 487 

[Iv2] V.I. Ivrii: Asymptotics of the ground state energy for atoms and molecules in 
the self-generated magnetic field, Preprint. larXiv:1112.5 538 

[IS] V.I. Ivrii and I.M. Sigal: Asymptotics of the ground state energies of large 

Coulomb systems, Ann. of Math. (2), 138, 243-335 (1993). 

[LLS] E. H. Lieb, M. Loss and J. P. Solovej: Stability of Matter in Magnetic Fields, 
Phys. Rev. Lett. 75, 985-989 (1995) 

[LS] E. H. Lieb and B. Simon: The Thomas-Fermi theory of atoms, molecules and 
solids, Adv. Math. 23, 22-116 (1977) 

[SW1] H. Siedentop and R. Weikard: On the leading energy correction for the statistical 
model of an atom: interacting case, Commun. Math. Phys. 112, 471-490 (1987) 

[SS] J. P. Solovej, W. Spitzer: A new coherent states approach to semiclassics which 
gives Scott's correction. Commun. Math. Phys. 241 (2003), no. 2-3, 383-420. 



33 



[SSS] J. P. Solovej, T.0. S0rensen, W. Spitzer: Relativistic Scott correction for atoms 
and molecules. Comm. Pure Appl. Math. Vol. LXIII. 39-118 (2010). 



34 



